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Abstract 

We have evaluated the contribution to the anomalous magnetic moment of the electron from six 
tenth-order Feynman diagrams which contain eighth-order vacuum-polarization function formed by 
two light-by-light scattering diagrams connected by three photons. The integrals are constructed 
by two different methods. In the first method the subtractive counter terms are used to deal with 
ultraviolet (UV) singularities together with the requirement of gauge-invariance. In the second 
method, the Ward-Takahashi identity is applied to the light-by-light scattering amplitudes to 
eliminate UV singularities. Numerical evaluation confirms that the two methods are consistent 
with each other within their numerical uncertainties. Combining the two results statistically and 
adding small contribution from the muons and/or tau leptons, we obtain 0.000 399 9 (18) (a/vr)^. 
We also evaluated the contribution to the muon g — 2 from the same set of diagrams and found 
-1.263 44 (14) (a/vr)^ 
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I. INTRODUCTION 



The anomalous magnetic moment {g~2) of the electron has played the central role in 
testing the validity of quantum electrodynamics (QED) since its experimental and theoretical 
discovery in 1940's l|, 

The precision of g—2 measurements has been improved steadily in subsequent sixty years 
, 4|. The Harvard group recently succeeded in measuring the g value of the electron with a 







substantially reduced uncertainty 
pubhshed in 2006 and in 2008 



jy using a cylindrical Penning trap. Their measurements 
6| are 



ae(HV06) = 1 159 652 180.85 (76) x 10"^^ 
ae(HV08) = 1 159 652 180.73 (28) x 10"^^ 



[0.66 ppb] , 
[0.24 ppb] . 



(1) 
(2) 



Taking the presence of the muon and tau lepton into account the QED contribution to 
the electron g — 2 can be written in the general form 



ae(QED) = Ai + A2{me/m^) + A2{me/mr) + A:>,{me/m^,me/mr) 



(3) 



where Ai can be expanded into power series in - 



A^4nf) + 4«(^)VA'"( 



a 



+ .... 



z = 1,2,3, 



(4) 



whose coefficients are finite calculable quantities, which is guaranteed by the renormal- 



izability of Q 




101, 
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D. Thus far the coefficients up to the eighth-order have been calculated 



due to hadrons 



21 



22 



23|, 



15 



16 



17 



18 



19 



20(1 . The small but non-negligible corrections 



2J| and weak interactions j25| are also known with sufficient 
precision. 

Combining the experiment and the theory, one can determine the value of the fine struc- 



ture constant a 



27[ 



a-\ae) = 137.035 999 084 (12)(37)(33) [0.37 ppb] , (5) 
where the uncertainties come from numerical errors in the eighth-order term 

HQ , an 

educated guess of the tenth-order term [28|, and the experiment ([2]), in that order. Note 
that, for the first time in three decades, the experimental uncertainty (0.33 x 10^'') has been 
reduced to a value smaller than the combined theoretical uncertainty (0.39 x 10^^). The 



uncertainty of this a is about 20 times smaller than those of other independent methods, such 



as a Rb recoil velocity determination in an optical lattice 
atom interferometry 



29j or a Cs recoil velocity in an 



3d . |31| . A new Cs measurement is now in progress, which is designed to 



obtain the value of a with the relative uncertainty 0.3 ppb [32j. Such forthcoming progress 
of the atomic physics experiments will enable us to check the validity of QED with the 
accuracy less than 0.1 ppb by examining consistency of various values of a. 

Turning back to the electron g — we find that the largest theoretical uncertainty now 
comes from the tenth-order term A^^^''. Clearly an actual value, not an estimate, of this 
term is urgently needed. There are 12672 Feynman diagrams contributing to Our 



on- 



37 



going effort to evaluate all of them has been reported in several articles 33, 



38 



39, 



34 



35 



36, 



401 ]. In this paper, we report the contribution from the diagrams belonging to the 



gauge-invariant set Set I(j). These diagrams contain the eighth-order vacuum-polarization 
diagram formed by two light-by- light scattering diagrams connected by three photons, which 
was constructed first time in this work. Although the Set I(j) consists of only six Feynman 
diagrams and it turns out to be numerically very small, it has features not found in other 
12666 diagrams contributing to the tenth-order electron g — 1. Thus it deserves a special 
treatment as is described in this paper. 

The primary purpose of this paper is to report the contribution of the gauge-invariant set 
Set I(j) to the mass-independent term ^[^^^ of the electron g — l. The contributions to J^2^^ 
from closed loops of electrons, muons and/or tau leptons are evaluated and reported sepa- 
rately in Sec. HVl Summing up all contributions, we obtained the tenth-order contribution 
from Set I(j) 



a(i°)(Set I(j)) = (AS^°)(Set I(j)) + Ve/m,)(Set I(j)) 



vr/ 



0.000 399 9 (18) ( - 



(6) 



The contribution from the tau lepton is smaller than the uncertainty quoted here. 

The contribution of Set I(j) to the muon g—1 can be obtained by replacing the external 
(or open) electron line by a muon line, keeping the internal fermion loops intact. The result 



3 



LLp 



LLc 



FIG. 1: Eighth-order vacuum-polarization diagrams LLp and LLc. There are two diagrams of LLp 
type and four diagrams of LLc type. When inserted into a photon hue of the second-order vertex 
diagram, they give the tenth-order diagrams contributing to the lepton g — 2 called Set I(j). 

of numerical integration gives the mass-dependent term of the muon g — 2 



The contribution from tau lepton is of order of the uncertainty quoted here. 

In Sec. HI] we describe how to construct the eighth-order vacuum-polarization function of 
Set I(j). Three possible ways are considered. In Sec. [111] two of three methods are described 
in detail. The utility of the Ward-Takahashi identity applied to a vacuum-polarization dia- 
gram and a light-by-light scattering diagram is particularly emphasized. Once the vacuum- 
polarization function is constructed, its contribution to the tenth-order anomaly is easily 
calculated. The details of the numerical results are presented in Sec. HV] Sec. |V] is devoted 
to conclusion and discussion. Appendix E] describes new features of the vacuum-polarization 
function for the diagrams of Set I(j) and also shows how to obtain its contribution to the 
magnetic moment which does not rely on the photon spectral function explicitly. An example 
of the structure of the integrand used in the Method C is shown in Appendix [B] 



a 



(i°)(Set I(j)) = (Af )(m^/me)(Set I(j)) + 4'°V/.Mr)(Set I(j)) 
+A^^^\m^/me,m^/mr){^ei I(j))j 




(7) 
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II. EIGHTH-ORDER VACUUM-POLARIZATION DIAGRAMS WHICH CON- 
SIST OF TWO LIGHT-BY-LIGHT-SCATTERING SUBDIAGRAMS CONNECTED 
BY THREE PHOTONS 



Two approaches are found in the hterature for deahng with the insertion of a vacuum- 
polarization diagram in the photon hne of the second-order vertex diagram. One is to take 
advantage of the fact that such an insertion amounts to replacing the photon line by a sum of 
massive vector particles weighted by the spectral function, which is the absorptive part of the 
vacuum-polarization function. Another is to insert the vacuum-polarization function itself 
obtained by the Feynman-Dyson rules. The first method is very convenient if the spectral 
function is known exactly [41], or in good approximation 42]. In most cases where such 
a spectral function is not available, however, one is forced to choose the second approach. 
The tenth-order diagrams of Set I(j), which consist of eighth-order vacuum-polarization 
functions inserted into the second-order vertex diagram, belong to the latter. This approach 
was initially developed in Refs. 43, 44]. (See also Eqs. (5.6) and (5.8) of Ref . [4^). A more 



transparent and compact form is presented in 461 ]: 



(8) 



where a*^^"*""^ stands for the (2+r;,)th-order electron anomaly that is obtained from the second- 
order vertex diagram in which the renormalized nth-order vacuum-polarization function 11*^"^ 
is inserted. A derivation of Eq. ([8]) is given in Appendix \M 

In the second approach the problem is thus reduced to an explicit construction of 11 
from the gauge- invariant set Set I(j) of Feynman diagrams. When twisted and flipped 
appropriately, two of the vacuum-polarization diagrams of Set I(j) (called LLp) are reduced 
to planar form with three uncrossed photons, and four of them (called LLc) have lower two 
of the photon lines crossed (see Fig. [1]) . Applying Feynman-Dyson rules formally to one of 
the LLp-tjpe diagrams we obtain 



KIM 



1 

(27r)8 W 
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where each closed lepton loop contributes a factor —l,pi are linear combinations of loop 
momenta h, I2, h, h and external momentum q, which enters at the n vertex and leaving 
at the u vertex (see Fig. [T]). 

The second LLp-tjpe diagram is obtained by reversing the direction of the arrow of lepton 
lines in the second trace of Eq. Q- By charge-conjugation invariance of QED it is equivalent 
to the first one. The LLc-type diagrams are obtained by exchanging 7/3 and 7^ in the second 
trace of Eq. iQ. All four diagrams of LLc-type are equivalent to each other. 

Of course formal expressions such as Eq. (Q are UV-divergent and meaningless until they 
are regularized properly. We follow the standard procedure to extract physical information 
from the expression and a similar one for LLc: 

(i) Make them convergent by the Pauli-Villars regularization of lepton loops and the Feyn- 
man cutoff of photon propagators. 

(ii) Renormalize them by subtractive renormalization, where subtraction integrals must be 
regularized in the same way as in (i). 

(iii) Remove the regularization terms from the final renormalized formula. 

These steps ensure that individual integrals obtained are finite. However they still contain 
terms which are not gauge- invariant. These terms cancel out only after they are summed 
over the gauge-invariant set of Feynman diagrams. Some details of the steps are described 
in the following. 

The integral (Q has eight UV-divergent subdiagrams, including itself. They are, namely, 
the hght-by-light-scattering subdiagram L (formed by a closed loop of lepton lines 1, 2, 3, 
4), another light-by-light-scattering subdiagram R (formed by a closed loop of lepton lines 
5, 6, 7, 8), a sixth-order vertex diagram V (formed by lepton lines 1, 2, 3, 4, 6, 7 and photon 
lines a, b, c), another sixth-order vertex diagram W (formed by lepton lines 2, 3, 5, 6, 7, 8 
and photon lines a, b, c), diagrams of the type [L in V) and {R in W), and the diagram 
itself that consists of all lepton lines 1, 2, 3, 4, 5, 6, 7, 8 and all photon lines a, b, c (see 
Fig. [2]). One more type of UV-divergence caused by L and R together generates no terms 
which contribute to the anomaly. Thus it can be ignored. 

UV divergences coming from L and R are only logarithmic and can be controlled by 
the Pauli-Villars regularization of the lepton loop. Control of UV divergences of V and W 

6 



requires Pauli-Villars regularization as well as Feynman cut-off of virtual photon momenta. 
In the latter, the photon propagator with momentum k is regularized as 



1 r^" 1 



dM'— — - , (10) 



P - A2 J ^2 (P - M2)2 

where the photon mass A and the UV cut-off A are introduced temporarily and to be put 
to zero and infinity, respectively, in the end. Finally, we must control the UV divergence 
involving all lepton lines and all photon lines. It is important to note that this divergence 
cannot be controlled by Pauli-Villars regularizations of two closed lepton loops alone. The 
quadratic behavior of this divergence comes mostly from three photons working together, a 
novel feature encountered for the first time in the eighth-order vacuum polarization. 

The sum n^'^(g) of all six diagrams, two of LLp type and four of LLc type, is gauge invari- 
ant and completely free of divergence after charge renormalization is carried out. However, in 
our numerical work which adopts the parametric integral formulation based on the topology 



of an individual Feynman diagram 



47 



48( 1 . it is more convenient to deal with the diagrams 



LLp and LLc separately. This means that we must go one step backwards and explicitly 
carry out the renormalization of logarithmic divergence from light-by-light scattering sub- 
diagrams, etc., as well as the quadratic divergence from the vacuum-polarization diagram 
as a whole. The logarithmic divergence is very mild and its removal by renormalization can 
be handled within the numerical framework keeping the gauge-invariance rigorously. 

The standard way to handle the quadratic UV divergence is to note that the Lorentz 
covariance dictates that Hg^, of either Q =LLp or LLc , consists of two scalar functions 

n^'^(g) = + q^q-'U^^^q') , (11) 

and note that Ug^lq) has the dimension of square of momentum so that the quadratic 
divergence (which is proportional to the cut-off momentum squared) is confined to the term 
proportional to g'^'^, or more precisely to the g-independent part of Il^g\q'^). Thus, if we 
write Ilg'\q'^) as 

n^fiq') = [n(f)(g2) - <)(o)] + n(f)(o), (12) 

the term within the parentheses is free from the quadratic divergence. 

As is well-known, gauge-invariance dictates that quadratic divergences in Ug^lq) of the in- 
dividual diagrams should disappear from the sum of the gauge-invariant set of the diagrams, 
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FIG. 2: Divergence structure of LLp. Subdiagrams are L, R, V, W, L in V, and R in W, and the 
whole diagram A of LLp. 
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and n'^'^(g) = 411'^^^ + 211^'^^ should satisfy the transversahty condition 

W{q) = {q^q'' - q^ g^"')U^^\q^) . (13) 

The scalar function 11*^^^ (g^) = 2n^''|p(g^) + 4n^'']^^(g^) defined by this equation is free from 
all subdiagram UV divergences. However it still has an overall UV divergence which must 
be removed by subtraction of n*^^''(0), which is nothing but charge renormalization. 

These observations lead us to three possible methods for obtaining the renormalized (not 
yet gauge invariant) amplitude ng(g^), which is the LLp or LLc part of 11'^^^ (g^) — 11^'') (0). 
The first method is 

Method A. Collect all terms of Il^'^{q) which are coefficients of q^q'^. 

Another approach is to note that Eq. f[T^ implies 

qx^'^iq) = 0, (14) 

which is valid for arbitrary q. Differentiating this equation with respect to q^ we obtain 

n'^'^(g) = -qx^U'-'iq), (15) 

in which one power of q is extracted explicitly. This has the effect of removing the quadratic 
UV divergence automatically. Thus, we can choose 

Method B. Collect coefficients of g'^g'^ or those of —g^^q^ of 11'^^ from the 
right-hand side of Eq. ( 1151) . 

Yet another approach is to start from the equation involving the second derivative of 
n^^(g): 

n^''(g) = JgAg.^^n^'^(g), (16) 

2 dq^ dqy 

which follows from Eq. (fTSjl and 

gAg.n^'^(g) = 0, (17) 

and symmetry of 11^^ in /i and v. Thus we may also start from the following rule in which 
two powers of g are extracted explicitly: 

Method C. Collect coefficients of g^g'^ or those of —g^^^q^ of fl'^'' from the 
right-hand side of Eq. ( 1161) . 

It turns out that Method C has a distinct advantage over the other two. Not only the 



quadratic divergence but also subdiagram logarithmic UV divergences, except for the one 
requiring charge renormalization, are eliminated as a consequence of the second derivative. 
Aside from this difference, however. Method A and Method B are equally useful and effective 
as Method C for carrying out numerical evaluation of the contribution of the Set I(j). 



III. CONSTRUCTION OF THE VACUUM-POLARIZATION FUNCTION U{q'^) 
A. Parametric representation of n(g^) 

Once the integral is made finite by regularization, we can safely deal with individual 
diagrams. In this article we adopt the method of parametric representation which has been 



successfully apphed to similar problems [45|. We begin by replacing the numerator, e.g., of 
Eq. ([9]), by an operator 

J^^" = Tr[7'^(]?i + m)7°(1^2 + m)7^(1^3 + m)7'^(i;Z)4 + m)] 

X Tr [YiV>5 + m)^aiV>6 + m)^f3iV>7 + ^hdV>s + m)], (18) 



where 



49| 



1 roo o 



and bring it in front of the momentum integration. (This may not be as straightforward as 
it sounds, and requires a more careful argument of Pauli-Villars regularization. But the end 
result is correct.) Then we combine all denominators with the help of Feynman parameters 
Zi, . . . , for leptons and Za, Zb, Zc for photons: 

pI pI pI fj- pi - 

= 10! [ (dz) « , (20) 

J {T^.=,^^{pl-^l)+Za{pl-X')+Zb{pl-X')+zM->^')y' 

where the photon mass A is introduced temporarily, to be put to zero in the end. 

As usual individual photon propagators may be regularized using the Feynman cutoff 
f llOp . Alternately one may regularize all three photons together as follows: 

1 

(Ell ^^iPl - mi) + zM - A^) + - A2) + zM - A2))ii 

= -11 / dM^ g ^ , (21) 

>^A2 (Ei=l ^iiPl -mi) + ZaPl + ZbPl + Z^pI - ZabcM'^y^ 

10 



where A is the UV cutoff and Zabc = Za + Zb + Zc. Let us assume that such a regularization 
is always done. Now we can carry out the momentum integration and obtain 

d% I d% [ d% I d% 



~ ~((ll-l)!)/((ll-9)!) J ^ ^^LPV^' ^^^^ 

where 

V = Zi234m^ + Z^QjsW? + Zabc>^^ - q^G, (23) 

^1234 = zi + Z2 + zs + Zi, etc., and 

G = -ziAi + ZaAa + 2:5^5, (24) 

assuming that the photon momentum q enters the diagram at the vertex /i, goes through 
lepton hne 1, photon hne a, lepton hne 5, and exits from the vertex v. Ai is the scalar 
current associated with the line i 45|]. Note that Ai is defined assuming that the arrow of 
fermion line 1 is opposite to the direction of q, whereas the photon line a and fermion line 5 
are in the same direction as q (see Fig. [1]). (Actually, any continuous path of q is equivalent 
to that of f l24l) . as far as the integral is made finite by regularization. Note that this may 
not be guaranteed for divergent integrals.) U is a Jacobian from the momentum space to 
the Feynman parameter space, (dz) stands for the eleven dimensional integration variables 
of Feynman parameters with the constraint that the sum of eleven Feynman parameters is 
unity. Although this integral is still logarithmically divergent, when it is regularized with 
the cut-off of Eq. fl2Tl) . the jF^'^-operation can be carried out correctly. (If necessary, we can 
introduce another cutoff parameter.) 

Collecting all numerical factors and bringing the operator JF^'^ back into the integral we 
obtain, for LLp, 



'^^vy^> - ^ (27r)8 Vtt/ ^'"•^((ll-l)!)/((ll-9)!)y U^V^ 

Before proceeding further we must carry out renormalization explicitly, following a well- 
estabhshed method. We mention here only few aspects that are specific to Set I(j). 

The first point to note is that the renormalization constants for the sixth-order vertices 
V and W are actually zero for the gauge-invariant quantity n(g^) by Furry's theorem since 

11 



there is no "self-energy" diagrams corresponding to V or W. However, they are nonvanishing 
for individual integrals, and must be subtracted explicitly from the unrenormalized integral. 

ogarithmic part of such a subtraction term can be readily obtained by the 



The leading 

i^-operation 45|, in which the UV divergent part of the standard on-shell renormalization 
constant is used. However, in the case of Set I(j), it turns out to be better to construct 
the exact and full on-shell renormalization term which enables us to avoid the trouble of 
calculating the residual renormalization term separately. 

Similarly, for the light-by-light-scattering amplitudes of L and R, we can avoid residual 
renormalization by defining the renormalization terms as the standard on-shell amplitudes 
defined with all its momenta external to it put to zero. The gauge invariant set of this light- 
by-light scattering amplitude is summed up to zero, which can also be shown by calculation 
with the dimensional regularization, hence no residual renormalization is needed. 

The integrals obtained by the Method A and the Method B can be shown to be analyti- 



cally identical using "Kirchhoff 's laws" on junctions and loops 45||. In the following we shall 
therefore consider only Method B and Method C. 

B. More on Method B 

We are now ready to consider Method B in detail. Let us write the integral for LLp 
symbolically, ignoring explicit multiple integration, as 

H = LSR, (26) 

where L and R are light-by-light-scattering diagrams introduced previously and 5* stands 
for the set of three photons connecting L and R. Then, the differentiation in Eq. ( fTSl) can 
be carried out as 

The first and third terms involve differentiation of the lepton propagators in the closed lepton 
loops while the second one is differentiation of the photon propagator. These differentiation 



can be carried out using the identities [45 1 



^ ^ -2I^^(^ + m) ,,^ , J (28) 



1 _ -2{p + qr _ ^29) 



dq/j. (p + qY {p + q) 
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This operation gives rise to an additional denominator factor which can be handled, for 
instance, as follows: 

^ n ' 



pI pI pI pI - f \ pi - mf 
(9111° 



n 



dml pI pI pI Ai pf - mf 
= 11! lidz) g . (30) 

J {Ell^^iPl-^l)+^a{pl-X')+Z,{pl-X')+zM-^W 

As a consequence, in Eq. ( l25l) . is replaced by —1/V^, 2! is replaced by 3!, —2D^' and 

—2{p + q)'^ are multiplied by Zi, etc., and then everything is multiplied by an overall factor 
—qx. Recall also that the direction of q is chosen to be opposite to that of p of the lepton 
line 1. In this manner we obtain 

= -| J {dz) (+2ziDf - 2.„D^^ - 2z,D^,) 

X IV [iijp^ + m)7"(i;32 + m)7^(I?3 + 111)7^(1^4 + m)] 

X Tr [Y{V>^ + m)7„(I?6 + 111)7^(1^7 + m)7c(1^8 + m)]:^. (31) 

Performing D-operation on l/V^ using Eq. f|T9l) . this integral can be expressed in terms of 
"building blocks" Zj, Bij, Ai, where i, j are indexes for lepton and photon lines. Of course 
it must be modified by various terms required for renormalization. 

C. More on Method C 

Let us now consider Method C. In this case it is more convenient to choose the graphic 
representation in which all three photon lines are parallel (or, uncrossed) in S. Then L 
can be replaced by a gauge-invariant sum of six light-by-light-scattering diagrams, which we 
denote as to indicate that it contains the vertex /j,. Similarly R is replaced by W^. The 
explicit form of of LLp is given by 

= 2[m"f'^{q, -p,, -p,) + U^^^'^iq, -p,, -pj + n'^^"^(g, -p„ -p,)], (32) 

where the light-by-light scattering tensor n'^"^^(g, —pa, —Pb, —Pc) is defined by 



n'^"^^?, -Pa, -Pb, -Pc) oc J dH Tr 



Y i 7° i 7^ i 7^ - 



^1 — m 1^2 — 111 — m — m 
pi = l -q, P2 = I - Pb - Pc, P3 = l - Pc, Pa = I, (33) 
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with the overall momentum conservation q = Pa + Pb +Pc- Actually this procedure gives six 
identical copies of the original six diagrams so that the result must be divided by 6. 

The differentiations in Eq. (1161) . where 11'^'^ is replaced by {1/6)L^SR'^ symbolically, can 
be carried out as follows: 

d d . d'L' dL' dS dL' dR^ 

dq^ dqy oq^oq^ oq^ oq^ oq^ oq^ 

dL^ ds ^„ . d^s .dSdR' 



dq^ dqy dq^dq^ dq^ dq 



dq^ dqu dq^ dq^ dq^dq^ 
Although this looks awful, it can be simplified greatly using Ward-Takahashi identities that 
hold for the gauge-invariant sum (or of light-by-light-scattering diagrams jsoj l: 

qxL^ = 0, q^R'' = 0. (35) 

Multiplying Eq. ( IMI) with qxq^j and applying Eq. ( l35l) . we obtain 

2 dqf, dqy 

The great advantage of this equation is that the derivatives like are UV-finite so 
that cut-offs can be safely removed and 11(5'^) can be evaluated without renormalization of 
subdiagram divergences. Of course the overall logarithmic UV divergence must be disposed 
by charge renormalization. 

Eq. fl36p provides the starting point of numerical evaluation by Method C. To perform 
numerical integration, one has to decompose it into non-gauge-invariant forms similar to 
LLp and LLc. This can be parametrized in the same manner as for Eq. fl3Tl) . in which 
light-by-light-scattering diagrams are treated as sub diagrams of the eighth-order vacuum- 
polarization diagram. 

In the following, however, we chose an alternate approach which emphasizes the gauge- 
invariant nature of the sets and R'^ of light-by-light-scattering subdiagrams. The set 

is a sum of six diagrams in which three photon lines and an external photon line n 
are attached to a directed lepton loop in every possible ways. We prepare another set R" 
similarly. We connect them by three photons to construct L^SR^ . This procedure can also 

14 



FIG. 3: Diagrams of L^SB!^ . Nine diagrams shown in the bottom three rows are obtained by 
inserting the external photons labelled ^ and v into the left and right fermion loops, respectively, 
of the diagram shown at the top. The (bold) internal photon line a carries the momentum q in all 
cases. Remaining 27 diagrams obtained by flipping the direction of fermion loops are not shown 
for simplicity. 

be stated as follows: suppose there is a diagram with two lepton loops connected by three 
parallel photons labelled by a, (3, Q from the top to the bottom, and we insert an external 
vertex /i into the left loop, another external vertex v into the right loop (see the top figure 
of Fig. [3]). Disregarding the directions of lepton loops for a moment, there are 9 ways of 
insertions. 

It is found that by flipping and twisting the diagrams three of them are topologically 
equivalent to LLp-tjpe diagram, whose photon lines a, b, c correspond to the cyclic per- 
mutations of a, P, (, namely, {/?,C;Ct}; and {(,a,P}, as shown in the second 
row of Fig. [3l Similarly, the remaining six diagrams are found to be equivalent to LLc-type, 
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whose photon hnes correspond to all six permutations, namely, {a, P, (}, {a, (, P}, {P, a, (}, 
{P, C, a}, {C, C(, P}, and {(, P, a}, as shown in the third and fourth rows of Fig. [31 

Next we consider the flow of external momentum q in the diagram. Three of four exter- 
nal photon momenta of a light-by- light scattering diagram are independent because of the 
momentum conservation. and in LLp are connected by three photons forming two 
loops in S. Two of three independent momenta of thus turn into two loop momenta 
and can be freely shifted. Therefore, to fix all external photon momenta of and those of 

in the L^SR^ we need to fix where the only one independent momentum q flows in the 
entire vacuum-polarization diagram. 

We define a fraction of the momentum q^ flowing in the line i, and introduce a coefficient 
di as 

qi = diq . (37) 

By momentum conservation, the sum of fractions flowing on three photon lines a, 6, and c 
of S must be equal to 1: 

da + db + dc = l . (38) 

Next we flx the flow of momentum through the gauge-invariant combination L^SR^ . We 
consider the following particular choices. 

Choice 1: q flows only through the photon line a. 

The combination L^SR'^ is decomposed into LLp-tjpe or LLc-tjpe diagrams listed in Fig. [3l 
in which q flows only on the photon labelled by a as shown by bold lines. For example, 
the diagram {a, P, (} corresponds to the LLp-tjpe diagram in which the fractions are given 
by da = 1, and df, = dc = 0, denoted symbolically as LLp {da = 1). Similarly, we can 
translate all nine diagrams of L^SR'^ variant to LLp or LLc with speciflc values of di. Then, 
the gauge-invariant vacuum-polarization function which contains two sets of gauge-invariant 
light-by-light-scattering subdiagrams is obtained by the combination 

^ [{LLp {da = 1) + LLp (4 = 1) + LLp (4 = 1)} 
o 

+ 2 {LLc {da = 1) + LLc {di, = 1) + LLc (4 = 1)}] , (39) 

where it is multiplied by 4 to account for the directions of lepton loops and divided by 6 to 
take account of duplicated copies. 
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Once we have selected the flow of q in the photon hnes, we can choose any flow on the 
fermion loops. For instance, we may choose the following flows for three LLp-tjpe diagrams: 

diagram with a speciflc g-flow fermion 1-4 fermion 5-8 

LLp {da = 1) di = -1, d^ = +1, 

LLp {db = 1) di = d2 = —1, d^ = de = +1, 
LLp {dc =1) di = +1, ds = —1, 

and other (ij's are zero. For each diagram, different routings of the external momentum q in 
the lepton lines is possible, but all give identical results. (This is nothing but a consequence 



of the "Kirchhoff's laws" for loops and junctions applied to LLj? j45l|.) 

In order to drop the unwanted terms, derivatives of S, second derivatives of L and so 
on, from Eq. (1341) . we must add up the contributions from three diagrams in the flrst line 
of Eq. ( l39l) . To do this, we deflne the derivative factor D'^'^ consisting of two /^-operators of 
Eq. ([28]) as 



three Lip-type diagrams \ j=l / \j=5 



A{-2ziz^D^^Dl - z^z^D^^Dl - z^z^Dt^Dl - z^z^D^^D^ - z^z^D^^Dl) . (40) 



This is derived by a consideration similar to the argument leading to Eq. (I3T1) from the 
structure of traces implicit in Eq. (136!) with the help of Eq. (!28l) . As seen in Fig. [3l LLp{dc = 
1) is in fact identical with LLp{da = 1), if the top and bottom of the flgure LLp{dc = 1) is 
reversed. Thus, we may double the contribution of LLp{da = 1) and drop LLp{dc = 1). Then 
the coefficient of ZiZ^D^D^ becomes —3 and no z^zgD'^D^ term is needed. The vacuum- 
polarization tensor of LLp-tjpe in Method C is thus given by 

X Tr [iiJpi + m)7"(1^2 + 111)7^(1^3 + m)7^(1^4 + m)] 

X Tr liiJp^ + m)7„(1^6 + m)7^(I?7 + m^dV^s + m)] . (41) 

where the indexes in D*^^''^^ are symmetrized with respect to /i and u. Eq. (14T|) is free from 
UV divergence except for the overall charge renormalization. The Pauli- Villas regularization 
is no longer required in Eq. (HTl) . 
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Similarly, we can construct 3'^" factor for LLc: 

B^" = A{-z^z^D^^Dl + ziz^DtDl + Z2Z^Di^Dl + z^z^D'^D^ + z^z^D>iDl) . (42) 

Another simple choice of the q flow is 
Choice 2: |g flows on all internal photon lines a, (3, and Q. 

In this case, all three of LLp- (or six of LLc-) types are indistinguishable. Thus, we find 
that the gauge-invariant set is 

^ {?>LLp {da = db = d,= 1/3) + 6LLc {da = db = dc = 1/3)} . (43) 

The D'^^ operators for LLp and LLc with this choice of q flow can be constructed in the 
same manner as those in Choice 1. The explicit forms of the integrand thus obtained are 
different between Choice 1 and Choice 2. Hereafter we shall call Method C with Choice 1 
and with Choice 2 as Method CI and Method C2, respectively. 

The /^-operators in Eq. fHTl) . etc., are applied to the \^-function on the right-hand side 



following the "contraction" rules 45| . Carrying out also the trace operations, the result can 



be written more explicitly in the form 

^LLp{q^ 



2\ f i'^^) 



TT 



(//(i),o + + {q')'H^i),2 + (g')'^(i),3 + {q^H^DA ^ 



+ {H(2),0 + g'^(2),l + (g')'^^{2),2 + {qYHi2),3' ^ 



-(i/(3),o + g'i^(3),l + (g')'^{3),2) ^ 



+(^(4),o + g^%),i)|yi^ 

H{5)fi . / A 



where the numerators Htr),h are expressed in terms of "building blocks" Bij, Zi, and Ai, i,j = 



1 



34 



45| and A is the UV cut-off in fl2T]) . A detailed structure of H(^r),h is presented 
in Appendix [Bl The charge renormalization can be trivially carried out, and the cut-off A 
can be put to infinity. We thus obtain the renormalized vacuum-polarization function 

nLLp{q^) = lim {nLLp{q^) - nLLp{0)) . (45) 

A^oo 

It is straightforward to translate Eq. (H4l) into numerical integration code in FORTRAN 



by carrying out algebraic manipulation involved with the help of FORM 
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5l|. 



IV. NUMERICAL RESULTS 



We are now ready to describe the numerical evaluation of the contributions of the Set I(j) 
to the electron g — 2. The largest contribution, which is mass independent, comes from the 
case where both fermion loops consist of electrons. We first made a preliminary evaluation of 



the coefficient of (a/vr)^ by VEGAS [5^ using relatively small samphng points. The results 
may be summarized as follows: 

AS^°^(Set I(j): Method B) = -0.072 8843 8 (57) + 0.0732 732 (118) 

= 0.000 388 9 (131), (46) 
AS^°^(Set I(j): Method CI g^"") = -0.059 537 (39) + 0.059 917 (30) 

= 0.000 380 (50), (47) 
AS^°^(Set l(j): Method CI gV) = -0.0161 450 (58) + 0.016 432 (161) 

= 0.000 287 (171), (48) 
AS^°^(Set I(j): Method C2) = -0.047 895 3 (60) + 0.048 290 3 (62) 

= 0.000 395 (87), (49) 

where the first and second terms on the first line of each case are from LLp and LLc 
respectively. Their vacuum-polarization functions are obtained from the terms proportional 
to g'^'^ in Method B, Method C2, and the first Method CI, and the terms proportional to 
qfj-qi^ in the second Method CI. Method B and Method C2 were evaluated on hp's Alpha 
station. For Method B, we used 10*^ points per iteration and 250 iterations, followed by 10^ 
points per iteration and 160 iterations for VEGAS integration. For Method C2, we used 
10® points per iteration and 100 iterations. Both of Method CI were carried out with 10^ 
points per iteration and 100 iteration on a PC with Intel's Core 2 processor. The Method B 
requires more sampling statistics for VEGAS than the Method C in order to reduce the 
uncertainty to the level of the latter, primarily because the renormalization is carried out by 
point-by-point cancellation between divergent pieces of the integrand. The individual terms 
of Method B, Method CI, and Method C2 are different from each other because of different 
treatments of renormalization and routing selection. Note also that the g'^'^ term and q'^q" 
term in Method CI give different integrands. The good agreement of these four cases within 
the numerical uncertainties provides a strong assurance of correctness of all calculations. 
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method 




contribution 


sampling points 
per iteration 


iteration 


B 




0.000 396 4 (59) 


108, lOio 


50, 100 


CI 




0.000 398 2 (31) 


109 


101 


CI 




0.000 393 8 (88) 


109 


50 


C2 




0.000 397 6 (25) 


109 


80 



TABLE I: A contribution to the mass-independent term of electron 5— 2, A\ , from the diagrams 
of Set I(j) calculated in various methods. All methods are analytically independent. The overall 
factor (^)^ is omitted for simplicity. The second column shows from which term, g^^ or q^q^ term, 
the vacuum-polarization function n(q'^) is obtained. The numeral in the parenthesis stands for the 
uncertainty in the last two digits. All calculations were carried out on RSCC. 



The most prominent feature of these calculations is that the LLp and LLc parts are nearly 
equal in magnitude and almost cancel each other. In view of similar analytic structures of 
these integrals, this suggests the possibility that cancellation takes place not only between 
the integrals as a whole but also between the integrands at many points in the domain 
of integration. If this is the case, one should be able to reduce the calculated uncertainty 
significantly by performing integration of the combination 2 x LLp + 4 x LLc. In order to 
verify this conjecture, we have carried out an extensive computation of the combination in 
both Method B and Method CI. 

The production job for evaluation of combined LLp and LLc were carried out on RIKEN's 
Super Combined Cluster System (RSCC) with 128 or 256 processors. It turns out that the 
vacuum-polarization functions n(g^) obtained from g^^j term and q^q'^ term have analytically 
different structures in Method CI and C2 even after combining LLp and LLc. This provides 
us with more opportunity to check our calculation. The numerical results obtained by 
various methods are summarized in Table [B 

Four values listed in Table [T] are independent of each other. Thus, combining these results 
statistically, we obtain 

A[^°^(Set I(j): combined) = 0.000 397 5 (18) (50) 

as the best estimate of the term y4^"'^°''(Set I(j)). 
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loop fermions 




contribution 


iteration 


(e Li) 

K^"! H'J 




2.281 (60) X 10"^ 


51 


( p III 




2 290 fll5) X 10~^ 


60 


( 11 ll) 

VP) H'J 




1.185 (17) X 10"'^ 


145 


( u. u) 

\H^ HJ 




1.284 (4) X 10"'' 


55 


(e 


y 


1.332 (95) X 10"^ 


100 


IP T 1 




1.62 (51) X 10"^ 


50 






4.988 (27) X 10-9 


50 


(/^> -t) 




5.007 (55) X 10^9 


50 






4.008 (99) X 10-1° 


90 


't) 




4.541 (13) X 10-1° 


50 



TABLE II: Mass-dependent contributions to the electron g — 2 from the diagrams of Set I(j) with 

/ \ 5 

e, /i and/or r lepton loops. All integrands are constructed in Method CI. The overall factor (^) 
is omitted for simplicity. The number of sampling points per iteration for VEGAS integration is 
10^ for all calculations. The numeral in the parenthesis stands for the uncertainty in the last two 
digits. All calculations were conducted on RSCC. 



The mass-dependent contributions to the electron g—2 involving muons and/or tau leptons 
are also calculated using the combined programs of Method CI with g^'^ term and with q'^q" 
term. Singular behavior of the integral caused by heavier leptons makes convergence of the 
integrand rather difficult. But, the contributions themselves are very small and currently 
of no interest compared with the experimental uncertainty. Therefore we do not need the 
precise values of the mass-dependent contribution. They are summarized in Table HTl 531] • 

Summing up all mass-dependent terms and the mass-independent contribution Eq. fl50l) . 
we find the total contribution to the electron g — 2 from Set I(j) given in Eq. ([6]). 

We also present the contributions to the muon g—2. They were calculated by replacing the 
external electron by a muon in the combined program of Method CI and/or C2. They are 
listed in Table IIIII The dominant contribution arises when both of light-by- light scattering 
loops consist of electrons. Statistically combining three results listed in Table IIIII of this 
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contribution, we found 



A^2^\m^/me){^et : combined) = -1.247 26 (12) 



'{e,e) 



(51) 



where the subscript (e, e) imphes that both fermion loops consist of electrons. Including 
all other contributions, we found the mass-dependent contribution to the muon g — 2 given 
in Eq. ([7]). By using the asymptotic expansion of the vacuum-polarization function with 
resp ect to the transfer momentum, Kataev obtained the leading terms of this contribution: 

54, m 



Af)(m^/me)(Set I(j) 



— a 



'2,1-1] 



+ 



{e,e) ■ asympt. 

5' 
' 4 



36 3^*^^'* 



+ 



--a 



[2,l~l] 



2.0237 + — 



TUp 



m, 



where 04^''"'' is the unknown constant term of the asymptotic expansion of the eighth-order 
vacuum-polarization function formed by two light-by- light subdiagrams. Comparing our 
result (ISTj) and the formula (!52l) . we obtain 



(52) 



[2,1-1] 



4.5529 + C ( — 



(53) 



V. DISCUSSION 

In this paper, we report the evaluation of the contribution to the electron g — 2 and 
muon g — 2 from Set I(j) which consists of six vacuum-polarization diagrams formed by two 
light-by-light scattering subdiagrams. The contribution to the electron g—2 given in Eq. ([6]) 
amounts 0.000 027 x 10^^^, which is far smaller than the current experimental uncertainty 
0.28 X 10^^^ in Eq. ([2]). Thus far we have no clear explanation of why the contribution from 
Set I(j) is so small compared with other diagrams of the tenth order. 

We have also demonstrated the utility of the Ward-Takahashi identity to deal with the 
light-by-light scattering subdiagram. Although it may not always help us to streamline the 
work-flow for writing numerical programs, we continue to examine its application to the 
computation of the tenth-order diagrams containing a light-by-light scattering subdiagram 
that have not been evaluated yet. 
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loop fermions 


method 


or Q'^g'' 


contribution 


iteration 


(e, e) 


CI 




-1.247 28 (25) 


50 


(e, e) 


CI 




-1.247 30 (15) 


71 


(e, e) 


C2 




-1.247 08 (31) 


25 


(e, At) 


C2 




-0.016 455 (71) 


20 


(e, r) 


C2 




0.109 9 (51) X 10-3 


20 


"t) 


C2 




0.149 (13) X 10-3 


20 


(r, r) 


C2 




0.189 3 (5) X lO-'' 


20 



TABLE III: Contributions to the mass-dependent term of muon g — 2 from the diagrams of Set 
I(j) with e, /X and/or r lepton loops. The overall factor (^) is omitted for simplicity. The number 
of sampling points per iteration for VEGAS integration is 10^ for all calculations. The numeral in 
the parenthesis stands for the uncertainty in the last two digits. All calculations were carried out 
on RSCC. 
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APPENDIX A: VACUUM-POLARIZATION INSERTION INTO THE 
ANOMALY INTEGRAL 

Analyticity of a vacuum-polarization function n(g^) ensures that once-subtracted disper- 
sion relation between its real part and imaginary part is given by 

Ren(q2) 1 p Imn(k2) 



1 f°° „^ Imn(k^) , 
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In the cases of vacuum-polarization of second, fourth, and sixth orders, the cut starts at 
g2 = 4m2, where m is the electron mass, and the imaginary part of n(g2) is nonvanishing 
for g2 > 477j2 only. For 11 (g^) of Set I(j), however, the cut starts at = because of three 
photon intermediate states, and is included in Eq. flAip . This is a novel feature encountered 
for the first time in the eighth-order vacuum polarization. Eq. (lAip also assumes that Il{q^) 
has no pole at = 0. This may be justified as follows: At the threshold = the 
absorptive part of n(g2) is proportional to the square of light-by-light-scattering amplitude, 
which is^ proportional to q^ because the light-by-light amplitude is known to be proportional 
to [56']. Meanwhile, three photon propagators could produce at most so that Imn(q^) 
behaves as or even higher positive power of g^ as g^ — > 0. Thus the singularity at g^ = 
cannot be a pole. 

Eq. flAip guarantees that 11 (g^) can be expressed by a spectral representation 



g2 Jq — g2 + A;2 
where 

- '-'-^ . (A3) 

The effect of inserting a vacuum-polarization diagram into a photon line with momentum g 
is thus obtained by replacing the photon propagator by a sum of massive vector propagators 
whose mass squared is k"^: 

This is easily translated into a Feynman-parametric integral formula. If a vacuum- 
polarization is inserted into a photon line Za-, we need to replace a photon mass by 
A;^, multiply the spectral function p{k'^) to the whole integrand, and integrate over a "pho- 
ton mass" k"^. This is the most efficient way to describe an effect of the vacuum-polarization 
insertion in the anomaly calculation. 

The spectral function, or the imaginary part of n(g2), however, is not always available, 
especially in higher-order cases 431]. On the other hand, we can directly construct n(g2) 
itself, or its real part, in the Feynman-parameter space using the FejTiman- Dyson rules. 
Thus our problem becomes how to express the effect of vacuum-polarization insertion by 
using the real part of 11 (g^). 
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Let us specifically consider the anomaly contribution from a diagram in which a vacuum- 
polarization diagram is inserted into the second-order vertex diagram. We will omit the 
overall factor a/vr for simplicity and set the electron mass m to unity. It is given as the 



integral over the Feynman parameters 



M,p= r dPpik^) [ ^4$-^, (A5) 







f/2 V + Zak"^' 

where a Feynman parameter assigned to the photon line is Za and that to the fermion line 



is zi. The explicit form of Fq, V,- ■ ■ , etc. are [45|] 



(dz) = dzidza6{l - Zi - Za), U = l, Ax = \-jj 

Fo = ziAiil -A,), V = zi- ziAi + ZaX^ . ( A6) 
Comparing Eq. (1A5I) to Eq. (1A2I) . we find 

= / ^^?(-n(,')) (A7) 



with 



q^ = --. (A8) 

Zn 



Substituting the explicit forms in Eq. ( ]A6[) and identifying zi = y and Za = I — y, one find 
that Eq. ([ATI) becomes Eq. ([8]). 

Higher-order diagrams contributing to the magnetic moment may have V"^ or higher 
powers of V in the denominators. We can easily extend the above method to such cases. 
Namely, the effect of vacuum-polarization insertion into a photon line Za is expressed by 
replacing the denominator 1/V"' for n > 1 according to a following rule: 

This rule is used in our forthcoming papers dealing with the insertion of vacuum-polarization 
loops in the magnetic moments of fourth, sixth, and eighth orders. 

APPENDIX B: STRUCTURE OF THE INTEGRAND OF METHOD C 

The integral Eq. (jUj) generated by Method C is very lengthy containing more than 30,000 
terms. The integrand, however, can be shortened by observing its structure carefully. When 
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the term proportional to g'^" is projected out from Eq. (HT]) . and Dj in D^"^ must be 
"contracted" with other -Dfc|45|. Knowing it, we can organize H(^r),h in the form 

Hir),k = Yl u I ( E d^^^^'^k] ( E ^^•^^■^^0 

l<k<l<8 l\i=l J \j=5 J 

4 8 

+ h E c^i^ri ^ Bij , (Bl) 



where T(^r),h and are expressed in terms of "building blocks" Bij, Zi, and Ai |34|. Then, 
the number of arithmetic operations is dramatically reduced and the computational time 
becomes less than one tenth of the program without the above artifice. 
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